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The Two-Higgs-Doublet Model is considered, in its CP-non-conserving version. It is shown quanti- 
tatively how vacuum stability and tree-level unitarity in the Higgs-Higgs-scattering sector constrain 
the parameter space of the model. In particular, at high values of tan f3, the model violates unitarity, 
unless some of the Higgs bosons are heavy. In the regime of large CP violation in the neutral-Higgs- 
t-quark sector, which requires tan/3 < 1, the Yukawa coupling parameter space (determined by the 
neutral-Higgs-sector rotation matrix) is reasonably unconstrained. On the other hand, the corre- 
sponding neutral-Higgs-6-quark sector allows for large CP violation at tan/3 3> 1. However, here 
the model is more constrained: Significant CP violation is correlated with a considerable splitting 
among the two heavier Higgs bosons. 



I. INTRODUCTION 

The Two-Higgs-Doublet Model (2HDM) has been pro- 
posed as an extension to the Standard Model (SM), in 
part because it provides an additional mechanism for CP 
violation [H, 0, [H, 0] • Various experimental observations 
impose non-trivial constraints on it. For example, the 
B — B oscillations p, [y, 0] and Z — > bb decay width 8| 
exclude low values of tan (3, whereas the B — ^ X s ~/ rate 9| 
excludes values of the charged-Higgs mass, M H ± , below 
approximately 300 GeV [10|. Also, the precise measure- 
ments at LEP of the so-called p parameter constrain the 
mass splitting in the Higgs sector, and force the masses to 
be not too far from the Z mass scale . While these in- 
dividual constraints are all well-known, we are not aware 
of any dedicated attempt to combine them, other than 
those of [H [3. 

There are also theoretical consistency conditions. In 
particular, for vacuum stability, the potential has to be 
positive for large values of the fields We shall fur- 

thermore require the Higgs-Higgs scatte ring amplitudes 
to satisfy perturbative unitarity (T5 , ITtI Il8j |. Taken to- 
gether, these constraints dramatically reduce the allowed 
parameter space of the model. 

The unitarity conditions are traditionall y p hrased in 
terms of upper bounds on the Higgs masses fig , [k^ . The 
present paper is devoted to a study of the vacuum stabil- 
ity (or positivity) and unitarity conditions. These limits 
will here be seen in conjunction with the CP-violating 
Yukawa couplings. We will study how the CP-violating 
couplings are constrained by the stability and unitarity 
constraints, for various Higgs mass scenarios. The com- 
bination with experimental constraints will be considered 
elsewhere. 
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In our parameterization of the model, we emphasize 
the masses and mixing angles. The latter are closely re- 
lated to the Yukawa couplings, and thus somewhat more 
"physical" than the parameters of the potential, to which 
they are clearly related. 



II. THE MODEL 

The present study is limited to the 2HDM (II), which 
is defined by having one Higgs doublet ($2) couple to 
the up- type quarks, and the other ($1) to the down-type 
quarks [2fJ. 

We take the 2HDM potential to be parametrized as: 
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(2.1) 
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Thus, the Z 2 symmetry will be respected by the quar- 
tic terms, and Flavour-Changing Neutral Currents are 
constrained [2l|. We shall refer to this model (without 
the Ag and A7 terms) as the 2HDM5. The more general 
model, with also A6 and A7 couplings, will be discussed 
elsewhere. 

We allow for CP violation, i.e., A5 and rrv\ 2 may be 
complex. All three neutral states will then mix, 



RM 2 R T = M 2 dia „ = diag(M 2 , M 2 , M 2 ), 



(2.2) 



where .M 2 is determined from second derivatives of the 
above potential. The 3x3 mixing matrix R governing 
the neutral sector will be parametrized in terms of the 
angles a\, a 2 and a 3 as in pi [22I]: 
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where c\ = cosai, s\ = sinai, etc., and 

7T 7T 7T 7T 7T , 

-2<"i<2> ~2 <a2 -2' °^" 3 ^2- ( } 

(In ref. [13], the angles are denoted as a = ai, on, = ot2, 
a c = a 3 .) For these angular ranges, we have C{ > 0, 
S3 > 0, whereas s\ and s 2 may be either positive or 
negative. We will use the terminology "general 2HDM" 
as a reminder that CP violation is allowed. 

With all three masses different, there are three limits 
of no CP violation, i.e., with two Higgs bosons that are 
CP even and one that is odd. The three limits are 15]: 



E x odd 
H 2 odd 
H 3 odd 



a 2 ~ ±7r/2, ai, a 3 arbitrary, 
a 2 = 0, a 3 = 7r/2, ai arbitrary, 
ai—OL-i — 0, oi\ arbitrary. (2-5) 



These limits of no CP-violation are indicated in Fig. [T] 
(For future reference, we display in Fig. Q] the full range 
— 7r/2 < a 3 < 7r/2, as is required for the general case of 
non-zero Xq and A7.) 
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FIG. 1: Limits of no CP-violation shown in the 02-03 plane, 
with the CP-odd Higgs boson identified. 

In the general CP-non-conserving case, the neutral sec- 
tor is conveniently described by these three mixing an- 
gles, together with two masses (Mi,M 2 ), tan/3 and the 
parameter /i 2 = Rem 2 2 /(2cos/3sin/3). In fact, since the 
Yukawa couplings are compactly expressed in terms of 
these rotation matrix elements, 
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Hjtt 
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■ 175 sm(3R j3 ], 
175 cos f3R j3 ], 



(2.6) 



the angles provide a rather physical way to parametrize 
the model. 

From Eq. (|2~2j) . it follows that 



E 

k 



RkiM^Rkj. 



(2.7) 



In the general CP-non-conserving case, both (Ai 2 )i 3 and 
(A4 2 ) 23 will be non-zero. In fact, they are related by 



(M 2 ) 13 =tan/3(X 2 ) 23 . 



(2.8) 



From these two equations, (|2.7p and (|2.8[) . we can deter- 
mine M 3 from Mi, M 2 , the angles a = (ai, a 2 , a 3 ) and 
tan/3 [22|: 
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MlR 13 {R 12 tan p-Rn) + M 2 R 23 (R 22 tan/3-i? 2 i 



(2.9) 



.R33CR31 - #32 tan/3) 
where we impose Mi < M 2 < M3. 

III. EXTRACTING THE A 



As discussed above, in the 2HDM 5 , with ImAs ^ 0, 
the two masses Mi and M 2 will together with a and 
tan /3 determine M 3 . Providing also M H ± and /i 2 , all the 
A's are determined as follows. Since the left-hand side 
of (|2.7|) can be expressed in terms of the parameters of 
the potential (see, for example, [Hj]), we can solve these 
equations and obtain the A's in terms of the rotation 
matrix, the neutral mass eigenvalues, fj, 2 and M H ± : 



Ai 
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[44 M i + ( Cl s 2 s 3 + s lC3 ) 2 M 2 



C 2 pV 2 



+ (cis 2 c 3 - sis 3 ) M 3 - s ii ], 
A? = -i[s 2 c 2 M 2 + (cic 3 - sis 2 s 3 ) 2 M 2 



S 2 pV 2 



\ 3 = 



+ (cis 3 + sis 2 c 3 ) 2 M 2 - cjfj, 2 
1 



(3.1) 



(3.2) 



r{ciSi[c^M 2 + (s 2 s 2 -^)M 2 ^ 



+ (44 - 4)M 2 3 ] + S2 c 3 s 3 (c 2 - st)(Ml - Ml)} 



1 



~2 [2M^± — fir] , 



(3.3) 
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A 4 = -\s\Mi + cMMi + 44M 2 +f i 2 - 2M 2 H± ], 

(3.4) 

Re A 5 - -i l-s 2 2 M 2 - 4sjM 2 - ^c\M 2 + M 2 ] , (3.5) 



ImAn = 



-1 



-{c /3 [CiC 2 S 2 M 1 2 - C 2 S 3 {C!S 2 S 3 + SiC 3 )M. 



+ c 2 c 3 (sis 3 - cis 2 c 3 )Mf] + sp[s!C 2 s 2 M 2 (3.6) 

+ C 2 S3(ciC3-SiS 2 S3)M|-C 2 C3(ciS3 + SiS 2 C 3 )M 3 2 ]}, 

where cp = cos/3, sp = sin/3. 

While Mf is given in terms of M 2 , M|, R and tan /3 by 
Eq. (|2.9p , it is more transparent not to substitute for M| 
in these expressions (|3.1|) (|3.6p . These equations are the 
analogues of those of [23j for the CP-conserving 2HDM5. 
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A. Large values of /i 2 

At large fi 2 > Mf, it is seen from (j3~Tj) and (|3"T2)) 
that Ai and A2 will eventually turn negative. This would 
violate stability and the model would break down. Thus, 
for fixed Higgs masses, there is an upper limit to /1 2 . 



B. Large values of tan/3 

According to the constraints of unitarity, reviewed in 
Sec. IIVC| the couplings Ai, A2 and | A3 1 cannot be too 
large. At large values of tan/3, where cp = cos (3 — > 0, 
the coefficients in Eq. (|3.ip multiplying M| and Mf will 
hence be constrained. When /1 2 is small, these coeffi- 
cients must be small. This requires |si| and | S2 1 both 
to be small. Otherwise, when /i 2 is relevant, the terms 
proportional to M| and Mf must balance against the 
yit 2 -term. 



IV. POSITIVITY AND UNITARITY 

We shall project the constraints of positivity and uni- 
tarity onto the tan (3-M H ± plane. Such a projection of 
information from a six-dimensional space onto a point in 
the tan (3-M H ± plane can be done in a variety of ways, 
all of which will lead to some loss of information. How- 
ever, we feel that this loss of detailed information can 
be compensated for by the "overview" obtained by the 
following procedure: 

1. Pick a set of neutral-Higgs-boson masses, (Ml, M 2 ) 
together with /i 2 . 

2. Scan an N = n\ x n 2 x 713 grid in the (ai,a 2 ,a3) 
space, and count the number j of these points that 
give a viable model. (Alternatively, one could scan 
over N random points in this space.) 

3. The ratio 

Q=j/N, 0<Q<1, (4.1) 

is then a figure of merit, a measure of "how allowed" 
the point is, in the tan (3-M H ± plane. If Q = 0, 
no sampled point in the a — (ai, CK2, 0:3) space is 
allowed. Similarly, if Q = 1, they are all allowed. 
An alternative measure 

Q+=j/N+, Q+>Q, (4.2) 

counts in the denominator only those points N + for 
which positivity is satisfied. 

Of course the 2HDM, if realized in nature, would only 
exist at one point in this parameter space. However, we 
think the above quantities Q and Q+ give meaningful 
measures of how "likely" different parameters are. 



A. Reference masses 

We shall impose the conditions of positivity and unitar- 
ity on the model, for the different "reference" mass sets 
given in Table H] (and variations around these). For each 
of these mass sets we scan the model properties in the 
a space. From these reference masses, some trends will 
emerge, allowing us to draw more general conclusions. 



Name 


Mi [GeV] 


M 2 [GeV] 


[GeVf 


"100-300" 
"150-300" 
"100-500" 
"150-500" 


100 
150 
100 
150 


300 
300 
500 
500 


[±(200)*] 
[±(200) 2 ] 
[±(200) 2 ] 
[±(200) 2 ] 



TABLE I: Reference masses. 



These masses are inspired by the indication from LEP 
that there is a relatively light Higgs boson [24|, here de- 
noted H\. The others, 7J 2 and H3, are then presumably 
more massive, and do not directly affect the LEP phe- 
nomenology. As an alternative, we shall also briefly con- 
sider the case of two light Higgs bosons, with the third 
one considerably more massive (see Sec. IIXI) . 

B. Stability 

Let us first explore the effect of imposing vacuum sta- 
bility, or positivity. The positivity conditions can be for- 
mulated as (for a general discussion, see Appendix A of 

My- 

Ai > 0, A 2 > 0, 

A 3 + min[0, A 4 - |A 5 |] > -y/Xi\ 2 . (4.3) 

Actually, we will use the notion "positivity" to include 
also the non-trivial conditions Mf > and Mi < M3 [see 
Eq. (|2.9p ]. We scan the a parameter space as discussed 
above, and show in Table |TT] the fraction of parameter 
points that satisfy "positivity" , as defined above. 



A 4 


-(200 GeV) 2 





(200 GeV) 2 


"100-300" 
"150-300" 
"100-500" 
"150-500" 


30.8-31.0% 
30.8-31.0% 
30.8-31.0% 
30.8-31.0% 


30.8-31.0% 
30.8-31.0% 
30.8-31.0% 
30.8-31.0% 


24.0-28.0% 

26.2- 31.0% 

27.3- 29.7% 
28.8-31.0% 



TABLE II: Percentage Q of points in a space for which posi- 
tivity is satisfied. 



We shall henceforth refer to the set of points in the a 
space where positivity is satisfied, as a+. The fraction Q 
of points in the a space for which positivity is satisfied, is 
around 30%. (The range given indicates the lowest and 
highest values found when scanning over 0.5 < tan/3 < 
50 and 200 GeV < M H ± < 700 GeV.) We note that an 
upper bound for this fraction is 50%. This comes about 
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from the fact that for a given value of (3 + at, in the 
2HDM 5 , only positive or only negative values of a 2 are 
allowed, not both Thus, with < a 3 < tt/2, the 

sign of a 2 will be given by that of 8 + a x : 

< 8 + ai < \k : a 2 a 3 > 0, 

-\n < 8 + a x < : a 2 a 3 < 0. (4.4) 

For small and negative values of fi 2 , "most" of the exclu- 
sion provided by the positivity constraint is already con- 
tained in the conditions M 2 > and M 3 > M 2 , without 
the explicit conditions (|4 .3[) on the A's. The conditions 
on the A's provide the additional exclusion at positive 
values of fj, 2 that is evident in Table UTl 

C. Unitarity 

Perturbative unitarity in the Higgs-Higgs sector im- 
poses upper bounds on the |Aj|. These relations have the 
structure 

-^r- 1 J2 a i X i + VQ(*i)\ < 1 ( 4 -5) 

where the coefficients are of 0(1) and Q(\i) is a 
quadratic expression [l6|, U3, Ell ■ These constraints can 
be expressed as upper bounds on the Higgs masses 
They are conveniently formulated in terms of the differ- 
ent weak isospin and hypercharge channels [li ]. 

When /j 2 = (or negative) , we see from Eq. (|3.1|) that 
Ai will become large when tan/3 1. Thus, unitar- 
ity will at some point be violated. This is illustrated in 
Fig. [21 where we show in yellow where at least one point 
in a+ satisfies positivity and unitarity. 

V. HIGGS-MEDIATED CP VIOLATION 

When the Yukawa couplings contain both a scalar and 
a pseudoscalar term, as in Eq. (|2.6p , the exchange of 
Higgs particles leads to CP violation via an amplitude, 
which for couplings to b and i-quarks is proportional to 

Y^^R^^-fdM 2 ) (5.1) 

and 

^P=I> 2 i*;3^/ t (Mj), (5.2) 
rr 1 sin 8 

respectively. The function f q (M 2 ) is in general some 
loop integral that depends on the Higgs mass Mj. This 
CP-violating effect is most important when the Higgs 
masses are not too close. Otherwise there are cancella- 
tions among different contributions, due to the orthogo- 
nality of the rotation matrix, 

3 3 



Also, since the functions f q (M 2 ) decrease for high values 
of M 2 , the effect tends to be larger when the lightest 
Higgs boson is reasonably light. 

Let us therefore focus on the couplings of the 
lightest Higgs boson, Hi. For maximal CP vio- 
lation in the H\bb coupling, R11R13 sin/3/ cos 2 8 — 
\ cos a\ sin 2a 2 sin [3/ cos 2 j3 must be large, requiring 

ai ~ 0, a 2 ~ ±7r/4, tan/3 > 1. (5.4) 

Similarly, for maximal CP violation in the H\tt coupling, 
R12R13 cos/3/ sin 2 /3 = \ sinai sin2a2 cos /3/ sin 2 Q must 
be large, or 

ai~±7r/2, a 2 ~±7r/4, tan/3 < 1. (5.5) 

When the two heavier Higgs bosons have a similar 
mass, M = M 2 — M 3 , the expression (|5.2p . for exam- 
ple, simplifies because of (|5.3[) : 

^ P ~ R 12 R 13 ^ [MM?) - ft(M 2 )} ■ (5.6) 
sm p 

Thus, also in this case is the coupling of the lightest Higgs 
boson of special importance. 

These two conditions (|5.4p and (|5.5p will be studied 
in the following. Common to both of them is the re- 
quirement that a 2 — ±7r/4. Vice versa, there is no CP 
violation mediated by the lightest Higgs exchange when 
a 2 ~ or when a 2 ~ ±7r/2. Also, we note that the 
conditions (15. 4p and (15. 5p do not refer to a 3 . 

In the case when two Higgs bosons are fairly light 
compared to the third one, by orthogonality, it will 
be the couplings of the heavy one that determine the 
amount of CP violation. For maximal CP violation in 
the H 3 bb coupling, |i?3ii? 33 | sin/3/ cos 2 8 — \ — c\s 2 c 3 + 
s\s 3 \c 2 c 3 sin/3/ cos 2 8 must then be large, requiring c 2 
and C3 to be non-zero. More explicitly, for small \ai \ one 
must have a 2 ~ ±7r/4 and a 3 ~ 0. At the other extreme, 
for |«i I ~ 7r/2, one must have a 2 ~ and a 3 ~ 7r/4. Sim- 
ilarly, for maximal CP violation in the H 3 ti coupling, 
-R32-R33I cos/3/sin 2 (3 = \c\s 3 + sis 2 c 3 \c 2 c 3 cos/3/ sin 2 (3 
must be large, also requiring c 2 and c 3 to be non-zero. 
More explicitly, for small |ai| one must have a 2 ~ and 
a 3 ~ 7r/4. At the other extreme, for |ai| ~ tt/2, one 
must have a 2 ~ ±7r/4 and a 3 ~ 0. 

We shall now proceed to study to what extent these 
regions of large CP-violation are allowed by the stability 
and unitarity constraints. 

VI. ALLOWED REGIONS FOR ^ = 

The unitarity constraint can have a rather dramatic ef- 
fect at "large" values of tan 8 and M H ± . While the gen- 
eral constraints on the charged- Higgs sector exclude low 
values of tan 8 and M H ± (see @, MM, S i, El US HI ) , 
the constraints of unitarity exclude high values of these 
same parameters. Only some region in the middle re- 
mains not excluded. For (Mi,M 2 ) = (100,300) GeV 
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FIG. 2: Percentage of points Q+ in the a+ space that satisfy 
unitarity. Four sets of (Mi, M2) values are considered, as indi- 
cated. All panels: (j 2 = 0. Yellow region: Q+ > (positivity 
and unitarity satisfied). The contours show Q+ = 0, 20%, 
40% and 60% (upper panels) and 0, 5%, 10% and 15% (lower 
panels) . The line at M H ± — 300 GeV indicates roughly what 
is excluded by the b — > S7 constraint [T^] . 



and fi = [see Fig. [2], unitarity excludes everything 
above tan/3 ~ 7 (for any value of M H ±), and above 
M H ± ~ 700 GeV (for any value of tan/3). 

For M 2 = 300 GeV, the percentage of points in a+ 
space for which unitarity is satisfied, reaches (at low tan (3 
and low M H ±) beyond 60%, whereas for M 2 = 500 GeV, 
it only reaches values close to 20%. 

The domains in which solutions exist (Q + > 0) depend 
on /i 2 : For negative values of [i 2 , the region typically 
shrinks to lower values of tan/3, for positive values of /x 2 
it extends to larger values of tan/3 (see next section). 
However, the maximum values of Q+ (at low values of 
tan/3), are little changed. 

As discussed in Sect. IIII B[ for large values of tan/3, 
the allowed solutions get constrained to a region of |«i| 
and |ck2 I both small. This is illustrated in Fig. [3l where 
we show regions of allowed a± and a 2 , for tan (3 < 5 
(lower part) and 5 < tan/3 < 10 (upper part). The 
masses considered are {Mi^Mz) = (100,300) GeV and 
(100,500) GeV. No allowed solutions were found for 
tan/3 > 7. As M 2 is increased from 300 GeV to 500 GeV, 
the allowed region shrinks. 

As discussed in Sec. [V] the Yukawa couplings of the 
lightest Higgs particle to b and t quarks is large for tan /3 
low or high, with |ai| ~ ir/2 or 0, respectively, and \a 2 \ — 
7r/4 in both cases. These regions will be referred to as 
"regions of major CP violation" and are indicated by 
boxes labeled "6" and "f ' in Fig. [31 We note that the 
boxes labeled "6" are empty, the model does not give 




FIG. 3: Allowed regions in the «i-a2 plane, for (M\,M2) = 
(100, 300) GeV and (100, 500) GeV, /i = and two slices in 
tan j3 as indicated. At higher values of tan {3 there are no al- 
lowed points (see also Fig. 0. Contours are shown at each 
negative power of 10, as appropriate. Yellow (light blue) in- 
dicates where the normalized distribution is higher than 10 -4 
(3 x 10~ 4 ); green (purple) levels above 10" 3 (3 x 10~ 3 ); red 
(dark blue) is above 10~ 2 (3 x 10~ 2 ). Regions of major CP 
violation are labeled "6" and "i". 



large CP violation in the bbH\ couplings for these mass 
and /i parameters. 

Since at high tan/3, |«i| and \a 2 \ are small, M3 will be 
almost degenerate with M 2 . The distribution of M3 val- 
ues is shown in Table [TTT1 It is seen that M3 is just barely 
larger than M 2 , in particular for high values of tan /3 and 
high values of M 2 (cf. M 2 = 500 GeV vs. 300 GeV). 



(M 1 ,M 2 ) = (1 00,300) (100,500) 




FIG. 4: Allowed regions in the 0:2-03 plane, for (M\,M2) = 
(100,300) GeV and (100,500) GeV, /i = and two slices in 
tan P as indicated. Contours and colour coding as in Fig. [3] 
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(Mi, M a ) = (100,300 [500]) GeV 



tan 13 


£<1.1 


1.1<(< 1.5 


1.5 < £ 


5-10 
< 5 


94.4 [98.5]% 
41.1 [82.5]% 


5.4 [ 1.5]% 
50.2 [17.5]% 


0.2 [ 0.0]% 
8.6 [ 0.0]% 



M,,M 2 , / u) = (1 00,300,200) (100,300,400) 



TABLE III: Distribution of M 3 values, £ 
and colour coding as in Fig. [3] 



Ms/M2- Contours 



is more 



The distribution in , of allowed solutions, 
spread out, as shown in Fig. |U 

At large values of tan /3 it turns out to be the isospin- 
zero, hypercharge-zero channel that is most constraining. 

When fi 2 < 0, the range in tan/3 is likewise limited, 
and the allowed regions in ai, «2 and CY3 are similar to 
those for the /i 2 = case. 



VII. ALLOWED REGIONS FOR Mi < fi 

The large-/x case is often referred to as the decoupling 
limit. It has received considerable attention in the CP- 
conserving case 27]. Within the framework set up by our 
choice of input parameters, it is natural to distinguish 
three mass scales: Mi, Mi and /i. Thus, there are three 
cases: 



{%) fi<Ah< M 2 , 

(a) Mi < [i < m 2 , 

(Hi) Mi < M 2 < /i, 



Sector! 

decoupling, 

decoupling. 



(7.1; 



If [i is "significantly" larger than Mi, the latter two both 
correspond to decoupling in the sense of Gunion and 
Haber [27j |, but from the point of view of CP violation, 
they can be rather different. 

In these regimes of Mi < fi, it is possible to keep 
Ai and | A.3 1 within the allowed range (not too large) by 
carefully tuning the other parameters. For fi suitably 
chosen (large), no part of the tan f3-M H ± plane is dis- 
allowed. From an inspection of Eq. (|3.ip for Ai, we see 
that |si| and/or | S2 1 must be small. But they can not 
both be zero, unless Mi is very close to p,. This region 
of small I s x I and/or |s 2 will also yield solutions for A3 
that are sufficiently small. In the following, we discuss 
the specific examples of (Mi, M 2 ) = (100, 300) GeV and 
(100,500) GeV, each of them for two values of /x. 



A. (Mi,M 2 ) = (100,300) GeV 

We display in Figs. \5\ and [H] the allowed regions in 
the a\-ai and a 2 -cv3 planes, for three slices in tan /3 
and two values of /1, "small" (200 GeV) and "moder- 
ate" (400 GeV). At tan/3 < 5 sizable regions in a space 
are populated with allowed solutions, but these regions 
shrink significantly for tan/3 > 5. 

From Figs. [5] and [H we see that for large tan/3 and 
increasing /x, the majority of solutions have values of a 2 




FIG. 5: Allowed regions in the a.\-a.i plane, for (Mi,M2) = 
(100, 300) GeV, n = 200 GeV and 400 GeV and three slices in 
tan/3 as indicated. Contours and colour coding as in Fig. [3] 
Regions of major CP violation are labeled "6" and "i". 



(M,,M 2 „u,) = ( 100,300,200) (100,300,400) 

7T/2I 





tt/2 





tt/2 



tar 




1 




tai 
5- 


r\ 


;< 


1 


tan£ 
<5 





-tt/2 Z\V2 -n 




FIG. 6: Allowed regions in the 02-03 plane, for (Mi, M2) = 
(100, 300) GeV, n = 200 GeV and 400 GeV and three slices in 
tan /3 as indicated. Contours and colour coding as in Fig. [3] 
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that move away from towards ir/2 (where CP is con- 
served, see Fig. [IJ. In this case of increasing fj,, M 3 will 
also increase, and the last two terms of (|3.1j) must com- 
pensate each other. In distinction from the case of /i = 0, 
\si\ and | S2 1 can then not both be small, one of them will 
approach unity, as seen from Fig. [5] According to (14. 3|) . 
A3 can not be too large and negative. From Eq. (|3.3p . 
this means that either si (i.e., a\) or s 2 (i.e., 0^2) must 
be large and positive, as seen in Figs. [5] and [6] 

There will for large values of tan/3 be a range of /it- 
values for which the allowed solutions accumulate around 
a 2 — ±7r/4. In these regions, the CP-violation in the 
iii&fe-coupling will be considerable. 

The distribution of M 3 -values is shown in Table IIVI 
Here, we note that for M\ < /i < M2, the values of M3 
are rather low (close to M 2 ), whereas for M 2 < /J. they 
tend to be considerably higher. 



(Mi, M 2 , fij = (100, 300, 200 [400 


) GeV 


tan/3 


€ < l.l 


1.1<(< 1.5 


1.5 < £ 


> 10 
5-10 
< 5 


58.6 [ 0.0]% 
41.8 [ 0.0]% 

29.7 [ 0.0]% 


40.0 [70.4]% 
55.3 [44.2]% 

56.1 [13.1]% 


1.4 [29.6]% 
3.0 [55.8]% 
14.2 [86.9]% 



TABLE IV: Distribution of M 3 values, £ = M 3 /M 2 . 



B. (Mi, Mi) = (100, 500) GeV 

For this case of "large" M2, we display in Fig. [7] the 
allowed regions in the ati— a 2 plane, for three slices in 
tan/3 and two values of fJL, 200 GeV and 600 GeV. 

In this region of large M 2 , we note that at high values 
of tan/3, the allowed regions in a.\ get constrained to 
values around ot\ ~ or ot\ ~ ±7r/2 with \a 2 \ increasing 
with /i from to 7r/2. (We recall that the a 2 — > tt/2 limit 
represents the case when the lightest Higgs particle, Hi, 
is odd, and there is no CP violation.) It follows from 
Eq. (|2.9[) that when |oi and \a 2 \ are both small, M3 is 
close to M2. Furthermore, it follows from (|3.1|) that in 
this limit, we have 



Ai 



2*3 JW 2 



2 c 3 m 3 



2 2i 
] 



(7.2) 



This should not get too large, in order not to spoil uni- 
tarity. Since M2 and M 3 are comparable in this case, the 
distribution in 0:3 becomes wide. This is analogous to 
the situation shown in Fig. [51 left panel. 



7T 

^2 



M,,M 2 ,/i)=(1 00,500,200 ) (100,500,600) 
tan/? 

>io U 

si 



-7T/ 



7T. 



-n/2 



tan/? 
5-10 tb 




eS 1 - 




-7T/2 7T/2 -7T/2 7T/2 



FIG. 7: Allowed regions in the ai—as plane, for (Mi, M2) = 
(100, 500) GeV, fi = 200 GeV and 400 GeV and three slices in 
tan/3 as indicated. Regions of major CP violation are labeled 
"6" and "t". 



As discussed above, when tan/3 ^> 1, the allowed range 
of M3 values get squeezed to a narrow band just above 
M 2 . This is illustrated in Table El 

When fi increases still, the situation becomes reminis- 
cent of that shown in Fig. [5l right part. The majority of 
allowed solutions move towards ai ~ and a 2 ~ ±7r/2, 
with small islands of additional solutions at ai ~ ±7r/2 
and a 2 ~ 0. 



(M 1 ,M 2 ,fJ.) = (100,500,200/[600]) GeV 



tan /3 


f <1.1 


1.1<(< 1.5 


1.5 < £ 


> 10 
5-10 
< 5 


93.7 [0.0]% 
84.9 [0.0]% 
74.1 [0.0]% 


6.3 [94.8]% 
15.1 [93.7]% 
25.9 [91.5]% 


0.0 [5.2]% 
0.0 [6.3]% 
0.0 [8.5]% 



TABLE V: Distribution of M 3 values, £ = M3/M2 



VIII. CP-CONSERVING LIMITS 

In addition to the general criteria for limits of no CP vi- 
olation given in Eq. (]2.5p and Fig. LT1 there are important 
limits in which there is no CP violation in the Yukawa 
couplings involving the lightest Higgs boson: bbHi and 
UHi. 
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A. CP-conserving bbH\ coupling 

The regions of CP-invariant bbH\ coupling require ei- 
ther Rn ~ (implying oil ~ ±7r/2 and/or a 2 — ±7r/2) 
or i?i3 ~ (implying a 2 — 0). These limits are shown 
in Fig. [3 We see from Figs. 02 M and [7] that both these 
categories of CP-conserving regions exist for fi = as 
well as for [i > 0. 

B. CP-conserving tii/i coupling 

The regions of CP-invariant ttH\ coupling require ei- 
ther R\2 ~ (implying a\ ~ and/or a 2 ~ ±7r/2) or 
i?i3 ~ (implying a 2 — 0). These limits are shown in 
Fig. [H We see from the lower panels in Figs. El [5] and [7] 
that such regions exist for [i — as well as for /i > 0. 



Limits of no H, CP violation 







-tt/2 



H, < 

bbH,, ttH, 


)dd 

X 
X) 

bbH,, ttH, 


X) 

Hi < 


X 

X) 
X) 

)dd 



-tt/2 ai tt/2 



FIG. 8: Limits of no CP-violation in the bbH\ and ttHi 
Yukawa couplings, shown in the ot\-ot2 plane. Also indicated, 
are the limits of Hi odd, where there is no CP violation in- 
volving any of the three Higgs bosons. 



IX. TWO LIGHT HIGGS BOSONS 

In this section, we report on some results obtained with 
Mi = 100 GeV and M 2 = 150 GeV (or 200 GeV). There 
are two questions: (1) Which parts of the t&n f3-M H ± 
plane are populated by allowed solutions, and (2) to what 
extent do such models provide CP violation? Concerning 
the latter question, we recall that two light Higgs bosons 
will to some extent act coherently, with Yukawa strength 
proportional (by orthogonality) to that of the heaviest 
one, H3. 



A. (i = 

This case is similar to the case discussed in Sec. IVIl in 
the sense that for /j, = (or small) , there is an upper limit 
to tan/?. That limit is lifted as fi is increased. At low 
tan/3 and /i = 0, for example, the majority of solutions 
fall in the domain of small |ai|, small |a 2 |, and large 013, 
as is required for significant CP violation in the i-quark 
sector (see the discussion at the end of Sec. [Vj) . At larger 
values of tan/3 (and fi = 0), the solutions still populate 
the small \ai\, small |a 2 | region, which is unfavorable for 
CP violation in the 6-quark sector. 



B. M 2 < n 

Let us first consider the case of tan/3 = 0(1). As 
compared with the case /i = 0, when fi increases, |ai|, 
I as I and \as\ all tend to move towards larger values, 
|d£i I — » 7r/2, |a 2 | — > 7r/2, as — > 7r/2. This limit does 
not satisfy the conditions for major CP violation in the 
i-quark sector. At larger values of tan j3, the solutions 
move towards intermediate and negative values of a±, in- 
termediate values of a 2 , and intermediate values of 0:3, 
which is a favorable parameter region for CP violation in 
the 6-quark sector. 



X. CONCLUDING REMARKS 

We have made a survey of parameter regions of large 
CP violation in the Two Higgs Doublet Model. Because 
of the many independent model parameters, it is difficult 
to extract a simple picture. For the admittedly limited 
set of parameters studied, it was found that considerable 
CP violation can easily occur in the ttHi coupling at low 
values of tan j3. In order to have significant CP viola- 
tion in the bbH\ coupling, on the other hand, it appears 
necessary to have a large mass splitting among the two 
heavier Higgs bosons, H 2 and H 3 . 
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